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The purpose of the research project ’’Multi-Disciplinary Optimization of Aeroservoelastic 
Systems” in its third year (October 1991 to September 1992) was to continue the development 
of new methods for efficient aeroservoelastic analysis and optimization. The main targets in 
this year were to complete the development of analytical tools for the investigation of flutter 
with large stiffness changes, to continue the work on efficient continuous gust response an 
sensitivity derivatives, and to advance the techniques of calculating dynamic loads with 
control and unsteady aerodynamic effects. 

An efficient and highly accurate mathematical model for time-domain analysis of flut- 
ter during which large structural changes occur was developed in cooperation with Caro 
D. Wieseman of NASA Langley. The model was based on the second-year work Modal 
Coordinates for Aeroelastic Analysis with Large Local Structural Variations recent y ac- 
cepted for publication in the Journal of Aircraft. The flutter time simulation work, Time 
Simulation of Flutter with Large Stiffness Changes”, was presented by Ms. Wieseman at 
the 33rd Structures, Structural Dynamics and Materials Conference in Dallas, Texas, Ap 
1992. The paper was also accepted for publication in the Journal of Aircraft. 

The work on continuous gust response, together with Arie Zole, a Technion Master stu- 
dent, has been completed. An abstract of the paper "Continuous Gust Response and Sen- 
sitivity Derivatives Using State-Space Models” was submitted for presentation in the 33rd 
Israel Annual Conference on Aviation and Astronautics, February 1993. The abstract is 
given in Appendix A. A full paper will be shortly submitted for Journal publication. The 
work extends the optimization model to deal with continuous gust objectives in a way that 
facilitates their inclusion in the efficient multi-disciplinary optimization scheme. 

Currently under development is a works designed to extend the analysis and optimization 
capabilities to loads and stress considerations. The work, together with Eyal Presente, a 
Technion Master student, is on aircraft dynamic loads in response to impulsive and non- 
impulsive excitation. The work extends the formulations of the mode-displacement and 
summation-of-forces methods to include modes with significant local distortions, and load 
modes. An abstract of the paper "Structural Dynamic Loads in Response to Impulsive 
Excitation”, submitted for presentation in the International Forum on Aeroelasticity and 
Structural Dynamics 1993, Strasbourg, France, May 1993, is given in Appendix B. 

Another work performed this year under the Grant was "Size-Reduction Techniques for 
the Determination of Efficient Aeroservoelastic Models” given in Appendix C. The work will 
be published as a Chapter in Academic Press Volume 54, “Advances in Control and Dynamic 

Systems”. 
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Appendix A 

CONTINUOUS GUST RESPONSE AND SENSITIVITY DERIVATIVES 

USING STATE-SPACE MODELS 


Arie Zole and Mordechay Karpel 
Technion-Israel Institute of Technology 
Haifa,Israel 


INTRODUCTION 

An atmospheric flight vehicle is exposed to air turbulence that causes time- varying aero- 
dynamic loads. These continuous-gust loads, amplified by the aeroelastic response o e 
vehicle, may result in critical structural design conditions. The response of the automa ic 
flight control system to structural vibrations may also play an important role. The gust re- 
sponse is analyzed by aeroservoelastic models which include unsteady aerodynamic, control 

and structural dynamics effects. , . . . , 

The gusts caused by air turbulence are defined in statistical terms by their power spectral 

density (PSD) functions and root-mean-square (RMS) values. Aviation regulations define 
the RMS gust velocity values for which the RMS responses over the entire structure shouldbe 

analyzed and substantiated. Classical frequency-domain methods are used to calculate e 

structural response by first calculating the response to sinusoidal gusts, and the associated 
response PSD at many frequency points. Numerical integration is then used to calculate the 

Modern aeroservoelastic modeling techniques, which are based on constant-coeffici , 
time-domain, first-order (state-space) formulation 2 , opened the way for a significantly more 
efficient way to analyze continuous gust response. Augmentation of the aeroelastic states 
by a gust filter 3 , in a way that results in a linear system excited by white noise, facilita e 
a direct solution for the RMS response. The main difficulty of the state-space modeling is 
in the requirement for approximating the unsteady aerodynamic force coefficient matrices 
by rational functions of the Laplace variable. The resulting model contains aerodynamic 
states which represent the time lag in the development of the aerodynamic forces. T e 
minimum-state method 4 results in a relatively small number of aerodynamic states, which 
yields efficient aeroservoelastic stability and gust response analyses. 

The inclusion of gust response criteria in a simultaneous structural and control design 
process may be beneficial to the aircraft performance and its structural integrity. The state- 
space formulation facilitates analytical sensitivity derivatives of gust response parameters 
with respect to various structural and control design variables. These derivatives allow the 
efficient inclusion of gust response in the cost function of automated integrated optimization 
schemes 6 along with other aeroservoelastic criteria such as flutter, control effectiveness an 
control stability margins. 

The purposes of this paper are to apply the minimum-state modeling method to state- 
space gust response modeling of realistic design cases, to model the gust filter m a way 
that avoids numerical difficulties, to develop the associated gust response equations, to de- 
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velop expressions for analytic sensitivity derivatives, and to demonstrate the efficiency and 
accuracy of the suggested modeling scheme in comparison with classical methods. 

KEY EQUATIONS 

Modal representation of the structural dynamics, rational approximation of the unsteady 
aerodynamic coefficients, state-space realization of the control transfer functions and the 
introduction of a modified gust filter yield the closed-loop state-space equation of motion 

{i} = [A]{x} + {#„,}«; 

where the state vector include generalized structural displacements and velocities, and aero- 
dynamic, control and gust states 

{x} T = [ £ i x a X c X g ] 

w represents a white-noise process and {I?u>} has a single non-zero term in the {x fl } row of 
Eq. (1). Structural response and external loads can be defined by the output equation 

{: V } = (CIW (2) 

Solution for the state covariance matrix [X] is obtained by solving the Lyapunov equation 

[A] [. X ] + [X\[A} t = -{B w }{B w } t (3) 


The mean square of a response parameter is 


r 2 v = [C,][X\[Cj 


(4) 


where [CU is a row in [C]. The computation of the sensitivity derivative of <r v with respect 
to a design variable p starts with the definition of £[A] and £[C7 y ] ( given in the paper for 
various structural and control parameters). Differentiation of Lyapunov equation (3) yields 
a Lyapunov equation for the derivatives of [X\. 




Differentiation of Eq. (4) yields 

d 


a /, = [C,}l p lX)[C,f + 2§- p lC,}[X][C,f 


(5) 


( 6 ) 


The same Schur decomposition of [A] can be UBed for all responses and their sensitivity 
derivatives. 


SELECTED RESULTS 

The modeling method and the solution process are demonstrated with a twelve-mode 
model of a transport aircraft at subsonic flight. The aircraft has one control surface and is 
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exposed to vertical continuous gust. The control system relates the control surface deflec- 
tion to wing-tip acceleration through a second-order control law and a third-order actuator. 
The generalized structural properties and the doublet-lattice unsteady aerodynamic data 
base were calculated using the MSC/NASTRAN finite-element program. Minimum-state 
aerodynamic approximations were performed using the MIST code 4 . The state-space model 
has 23 structural states (the heave displacement state is eliminated to avoid a zero root), 

4 aerodynamic states, 5 control states and 3 Dryden gust states (with a low-pass filter to 
avoid a noise term in the output equation). 

Open-loop flutter and gust-response results are within 2% of those obtained by frequency- 
domain methods. The same aerodynamic data base was used for the state-space an 
frequency-domain solutions. The aerodynamic approximation (not required in the frequency- 
domain solution) took about 7 cpu seconds on a VAX 9000 machine. One approximation is 
used for numerous flutter, gust response and sensitivity derivative analyses. State-space cpu 
times for flutter and a single gust response solution were 0.44 and 0.12 seconds respectively. 
The cpu time for the respective frequency-domain solutions were about 15 times larger. 

The variations of wing-tip acceleration and wing-root bending moment RMS responses 
with the thickness of an additional plate near the wing root are shown in Figure 1. The 
expected values are based on sensitivity derivatives. It is clear that the the sensitivity 
derivatives are accurate. The additional computation time required to calculate a derivative 
is less than that required to calculate the response itself because both computations are 

performed with the same matrix decomposition. 

It can be concluded that the proposed modeling and solution procedures are suitable for 
inclusion in efficient automated design procedures with numerous response parameters and 
design variables (structural and control). 
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Appendix B 

Structural Dynamic Loads in Response to Impulsive Excitation 

Mordechay Karpel and Eyal Presente 
Technion - Israel Institute of Technology 
Haifa, Israel 


The dynamic response of aerospace structures to impulsive excitation such as store ejec- 
tion and hard landing often yield critical design load cases. Time simulation of the dynamic 
response is usually based on the modal approach where a subset of low-frequency normal 
modes serve as generalized coordinates. The dynamic load distributions may be calculated 
by either the mode-displacement (MD) method, which is based on the generalized displace- 
ments, or the summation-of-forces (SOF) method, which sums the inertial, aerodynamic 
and excitation forces. While MD loads are easier to formulate and analyze , SOF loads are 
usually more reliable and require a smaller number of modes. The proposed paper presents 
a comparative investigation of the MD and SOF loads due to local impulsive excitation, arid 
a new method which improves the accuracy and efficiency of the MD analysis dramatically. 
The new method is based on application of large fictitious masses at the excitation points 
when the modes are calculated. The fictitious masses cause local distortions which are re- 
tained in the modal information when the masses are removed in the response analysis. T e 
local distortions allow local response effects, which are shown to be essential for accurate 
loads. Fictitious masses are also used to generate artificial load modes which yield simp e an 
efficient expressions for integrated shear forces, bending moments and torsion moments at 
various wing sections. The time-domain models include unsteady aerodynamic effects which 
are based on the minimum-state method for rational approximation of the aerodynamic force 
coefficient matrices. The time simulation is based on analytic expressions. 

The methods are demonstrated for store ejection loads on a high aspect ratio wing. 
Various response cases with and without aerodynamics demonstrate the MD and SOF loads 
versus number of modes taken into account. Typical 10-mode MD integrated section-load 
errors without fictitious masses are 10% at wing sections far from the excitation point and 
30% at sections close to the excitation point. SOF errors in these cases are about 0.25 of the 
MD errors. With a single fictitious mass of about one third of the total wing mass loading 
the excitation point, the 10-mode MD errors are reduced to less than 3% at all the wmg 
sections. SOF errors with fictitious masses are similar to the fictitious-mass MD errors. The 
new method allows high-accuracy dynamic load computations with the simple MD method 
which does not require the complicated evaluation of unsteady 
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where the i indices relate to the ;th rows of the respective matrices and ^ u o‘ ~ 

where the number of unknowns in each solution is 3 (n m + n c + n ff ) -f n^. 
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The statically residualised model of Eqs. (47) and (48) can now be con- 
structed using Eqs. (40) and (51). It can be shown from Eqs. (10), (49), 
and (61) that the matrix coefficients of the resulting state-space equation 





0 ~ — o 


D • f * 

K g.»tsj * 

x W * 

«fl DO 5 '** 

V e « 


"3 *T3 C >» 

5 v d *o h: 

; s « 1 c E 

L ,£5 g T3 ^ 

‘X no 3 a 5 
r '*■* *72 _ o <w - 
^ D. o ' 

; jj 'S *5 S w 

,is2^ = 

0 o . S - - 

> c o _ 2 

! ^ x. m 

* O ** lO V _ 

" 3 JS _ % 

„ C 3 o ^3 E 
2 O D* cr 3 V o 
; aj c5 * ^ o 

ai .*_» - _ ^ U 

r.!2s » 1 

> 1 5 -'SI 3 
is 3 o - H t 

1 s~r s 

5 o *Cj O s 3 »o 

1 *J2 <— r' . 3 *fl rl, 

j 5 -o <s o* . 

* S c « -a v cr 
, £> « « -3 ^ w 

o s - 

0 u'Q * a e 

s E § 11 1 2 1 

* ,2 a.5 ff 

< E s d, 

2 C 3 w- 3 — o 

a 2 JS 60 S 2-2 

3 ~ .5 83 s T3 
r w v J> 

J? M JO 'v *2 V O 

1 r^S • S.2 -2 

£ J 1 a * * § « 
2'^ oil * js 
C ^ 33 <c d — 


^ »Q Jl 

7=r — <« 

% o ■£ 
o, i *5 
i S. I 


ii ii 


!*g g> % 3 

| •£ | £ 3 

JK •§ a g 
> to £ g S 

• c t 

E o 2 « £ 

So 5 .a £L 


=< £ •= 2 ^ -2 I 

2 2 s •! 5 a I 

^ ^ *sb s £ 55 « 

vn C Ii bj m 

0 u .2 s -® £ S 

TS -| >, g S u 

1 I -5 I ?5- 
6 & §•«« = § 
”3 *5 2 5 -= = -a 


S 2 ° 5 c -a 
.u 2 So—.- o g 

-- « _ oo ^ ^ 3 

E^-S-s §t 

JS « 1 .1 1 1 s 

“ W M c> 

« £ v .2 * • js 

£ -3 js g* id .S 

*£ H ° 5 J 

•5 ^ o ^Q o ' 

<— • w C w 

j:j as ^ ts e — T5 
-- 3 .2 £ -T3 -T3 

>i ; s'® £ » 

.2 3 g. O J3 g, ® 
os w o « -** 


Si m r* ~ 

o c • > £ « v 

•« 2 r^r — O m o 

-3 - £. — -5 s ► 

5 w •> 2 sr <3 
a E > g|io 

6 5 — b S *« S 

w T3 ^ g g; * c 

u o tli 5 2* iS £ 

•s sS ! ‘ H.2 


£ o co ■ u S OS 

i < < js *1 a 5 

s § .s s 

? «a *c 


^ • xe J5 

r» jj 1 — « P •> 

£ c M "0^-0 

H I i .» 1 1 a 

Q ^ 3 H ■** 



y a j, t; -o b c 

C cl £-'-;'** w 

oj P 3*rcs 5 u l — ■ 4>'~S 

e^-S^^gjeoMO 

-X3 — ” _ *« o ^ ^ ^ t 

s 5 5.^ o E 

* _ 3 £■£':: ” w u J" 

-- C) cr •-« 0 2 ci d p 
^ «nrv,^C^’ 3 ^»-«2 

5 11 i ^ 8-i;i 

►■•J-S.SgS* 55 
- is v 5 -g ; e - s *S 

OS r-C *7* r* _ , -, * M 0} 

: -*-» O C_4 ^ C O tc» U* 

j cd Cm ^ t ■ r T '} •-< m 

! o 5 » "2 -o 7 w .s "m 

; u ^sSs!” - •» »> 

. o „■ E.2 ? iii s “ 


o « E | § I 2 «- 1 | 

■«|S a 8.Y-g St* 

v d o lo i, ii o 1 - 
■■^ k o ■** - S^^rO 

P S. -e *b o ® t*, > o c; 


£• to c , — „ „ O 

‘■S’ § |& !l 5 E g-2 2 

> 3 “ * 7 * o^~ 

:gj='T3l^ES^^ c 

> -a cfl £ 65 tj £ j^ ~ ~o E 
og«i— ^g,-oc^ 

j^=.^tE w.^ c3 -n 


2 ►* t &■ J! h j 3 

J j: T3 T *, E S 

> W o si R 


: Jr ^ e 

'TIL 


T -d"g c “ 

£ -3 8 51 -is 8 gj 

to o S^OojO^c 


to o 
1 c « -i* ® 

‘ .3 41 LO 


-W qj 

ill «2 sl^ : 

o c *> „b Cu cj X 

b js to'Tj^i-, a-r ^ E 

ro*C®Eo • Ji I w: 

2 > « : -j « ; >1 o o 

2 -a oT^ g 5 © u m - 

' k g® s«r>is rt 

r- d a 'l % a *> 

. H •> 11 « P L ?. a -3 5 .n 



O 



.“ IS to ” £ o ^ 

E .2 o 5 £ >- 5 

m « tor; ® 

“ s LO o H -E u 

r 0 g =3 

■5°°:^ 3 13 
2 - -r 5 o c “ 

«s j* 0 > ~a J: <u 

s D 2 £ * c ~a 

S 5 - *3 ^ ^ a> 

5 5°J/ & 

C ^3 ^ Oo - 2 

O « r-> m 0 7 — a> 

Z. t;' Hj 

41 C> 4 > *0 

2«o 

J “ o £ 11 -3 « *3 

01 05 3 -* <y 

43 X r — co w 2 

td ^ «•>. Cfl M 


2 j-i B g w J s 
«> = 0 n \ .2 0 3 
II =3 m ^ ^ ^ 

o?!! i s=*’ o5 S 

*7 _ *j id w 'T on _ 

s!d= 5 2 s - 2 

»B rt 2bo8o/< 

-§ 5 ^ -9 2 ff ® ^ ^ 

S *? 2 

.£ 5 65 L -s 5 c ^ 

IboJsuE 16 ®^ 

£ • TS £ £. ^ - 0 
‘ c l= 0 -S’ " 17 

» H 8-2 2*3 §E 

®2|h^.^ s •= 
g .= o “ - 5 • S 

£-2|°^ £ J5 sp 
2 S ^ 1 o ? .2 s. 

P to t-n - *» cn ® 2 

t2 - *> to to sj 

«J C * c <0 «J -n 
-O O O C — • — ” 


A^» "o w ^ O « 

W Cl ~ — 5 


§ o J ^ P 

3 2 41 

C o u 2 -* ' 

5 .»: s. 
“3 H g-S 

5 J rt >> “ 

■*2 — - H C* 

c o ^ . o 

„ >• » r— -V , - 

‘X " 3 O"^ 

° «r> ti ^ O 

>< ‘B _o v -* 

•J = ^ 5 0 ' 

g = g s E 

E ^ 8 2 £ 
•3 o £ Z ^ 1 

c ° — .H -S 

0 . x c - . 

bO 9 t> 

(U = 3 B C 

•- 4) ,o C E 

.X *— f ‘w V 

«£ 2^5 

j; « £ J> •= 

O ^ jS *» 

T 0 2 .£ | 

J. g»J si 

01 rv * — » 7? 

*-■ -3 ^ O 

« E a - 0 
-r <5 o -j c 

^ *0 W ? — 


(s)/Cictn3s®I 






>> oJ 
£ T3 ^ 
o 

.5 S -- 

cj OJ 

dog 

5 £ S 
XI J ~ 

fll JJ - 

> CT3 

u, * V 

<n 3 

w £S o 
3 

O M 
W d <u 
D :g -Q 


*0 3 *w 

JS v 
^ 3 

Jc “ 2 

2 ’ 5 
JS -j 

-*■» ^ 4» 

UQ J3 


** O 

d . 

— ' 4, 


up 

** J -o 

>» 4> 

_ o 

o ^ 

XI JS 

d 


1- 09 

> 5 c 
« ® 3 

3 

U s 
£ 2 

O MB 

d jj 


- g -o 
• 2 .5 *> 
° E •- 

m .*2 m 

M 0) 

V ^ IM 

d S Ou 

1«7 

w s c 

V 


* t> . 

** v ts 

-b» JS <fl 

3 o, 0 

« E -a 

c d g 
JS K £ 
+* V 3 


d 

>> j 

-o 4> 5 


> * 5 

« w d 
-r3 

a o -6 

s £ s 
:js s 

•! •§ cr 

-2 d m 

22* « 

4) O 
T3 3 IS 
4> ~3 

■S t ."2 
S to - ® 

c d t>o 
t- ^ c 
CX .*« 

*> s -a 

•« .£> 1 ? 


5 ^-g 
« S £ 


£ -3 S 

= *o ° 

aj 4 J u, 


« 3 ,2P 

4> _□ 

-c d 

O — 4> 


3 -5 d 
cr ^ j= 

4> JS 


i .5 &, 

.2 2 5 


r ^ no 
o ^ u -x? 

lO O U 

II p 3 > 

- ° v -S 

^ « 

- JS *- g- 

* H 2 S 

m l- 

«J .3 

OO C ^ 

! £ © <u Q 

; d _ 3 

1 C ^ i T3 

~ c S c 

; q « « si 

1 ^ cs 

I <1 o „ *-3 

^ V ^ "N 


-*-* <6 
.5 a ^ 


B 2 a 

d rt ~c 


w .3 c . 

> ■** 4J 

— t> 3 _ 

? g Z 

• £ c*T 

ill* 

v 5 v 

r2 8 5 

8 JS S 

o «* . 


^ co" d 

- o >* 

5 w 

° cs" *5 

(O * 

.a S 

« to 3 

c »--< d 

S o s 

d --° 


« -fi 

►-a 

O js 

to 


j'oS 


« c •— 

X ° -S 
5 £s 

3 C 

d 5u 

c> cr 

5 o: ^ 
c /5 ^ g 
bo o 
5H E 


d o T3 0 m 
£ C . O 

5* d v 52 w 

CX c o> 

d - oo D 

5 « O t- 

■2 c o > 

5 *5 H * 
i . t-" *3 

-G to no O 
- G v o ^ 

i o -3 ° d 

^ o d - > 

6 "d g 5 to 

£ 1 ! S>°“ 

§ ? s°.S 

° S * g -a 

o d ^ R c 

H "O cm o d 


•s i 

o o> +» 
u a> a 

«> o U 
* ■ 
pG o tn 

-* „ G 

»-« .2 o 
° , w 

42 VT3 
3 D 
“ to 
£ ^ 3 

a a » 
« g S 

>> e « 

o v d 

c a*; 

4> Q d 
3 g 

cr •> 

Z> J3 


y ^ 

nO 3 

ri js c "5 
£ u d Z 

S “ co-1 

«§ B S g 

flog 1 ' 

„**■»> 

a°5 

5 s II * 


^ m 

^ fid 
^ o .t» 

£ ° * 
-D n 

«■- 5 C 

> *> O 

0 i 

”S .H -ST 

» E 5 

d d o 
to G 
m >* ^ 

^ £ - 

1 SQ 

d » d 


d *— so 

HS f § 

5 “ S u 

(4 5 JS 
3 — *5 o 

^ 4) ^ 

^ js w 2 

*5>- S £ 

•2 « S i 

.< < 

k«, — » c- 

O 

u m 3 C 

t * 5 5 
*> *5 « d 

Ji “ ^ 4> 

— w .5 -« 

3 4> i 

CS s ^ L 

» S 3 Tg 

.£ 0 3 .£ 

^ £ S * • 
2 M S 
S » M ° 

•> « -2 ■= 

- ** Q o 

w- JS O o 
O r- '* J 

S 4» 

C SS 3 
O " d -O 

•£*;*£ 
d — M 

2 , -S >. 2 

E £ r: 

O 3 3 V 
u o « 

^ .2 2 . 55 

< C cr - 

d £ ^ 

.o g-*T 3 

n E s - 
jf 8 S. < 


3 

to 

to OO pH O 
CO t— o O 
t— < o o o 

O" 

to 

0.63 

0.56 

0.35 

0.00 

3 

to 

2.60 
1.94 1 
0.00 
-0.27 

&■ 

to 

0 O CT> O to 

tT 0> O tO 
r-I O O O 

i i » 

3 

* to 

0.16 

0.22 

0.02 

-0.06 

tF> 

* W 

0.28 

0.31 

0.03 

-0.18 


0.8 

0.8 

0.0896 

0.0544 

Data 

»— CM CM CM 






g W £> S .= 5: 3 g '3 S' S 

pL rt bb 03 O /J 

a .i«ss; 3 «i-is 

S 2 •= - o .2 .S S g g * 1 
^ & I 11 5 1 3 a 3 S --=3 


M K ^ jo .a .s 

■§ £| “ iga 

C rt £ ^ p 


*- 5 1 X 2 1 §•“ a a = - 

~CJ jS 01 .72 -» O « rJ g^* >* .2 


.2 01 .72 o « rf R >? .2 
£ <0 ® CX o> fc- — h- 

O.E ^ ^ ® ri 72 O r; ft, 4) 

<-X Z* «co £ ^ ■*» c - c ^ ^ 

* «_* g ^ 2p l _o c;> ri*n 


S.'o.S .S5 u 0 °-fl as 

« " Sh # ?i; - S S 8 S 

SOc3 _£.Sgrt oE=3 

£ *3 «*J tc 3 "3 r- t- •- » 0) 


. -2 -3 E? £?-g 

5 5 « S 2 1 


.«> £ J! * t • 33 4? & 


?--if 5-| £^ I'-g &&-g 

-|o«> ? -ss^|||| 

isl^.in-^5 S S | 

^ Lj — CT' tr ^ > o> cj nl 

u- “ 3 ?L 72 t£ ^ > JS cJ 

0 * c •*■*• o zi ^ ^ c m 

03 ,- W) 4) r- 4) o n ft, .2 

1 J S .2 | »|«g “S^-f S 

"rt “3 --3 w .2 g £ § w JE 

d^‘£.2tfij2«) 3Poq | Jri 

■3 — ? c ® ^ > .s ^ w 


H 


S" 

03 

o 

,03 

. 

c 


. 

fc= 

(73 



bp 

d 

c 

o 

03 

_o 

*S 

03 

*wu3 

^-a 

4J 

ed 

tz 


> 

RJ 

01 


15 
1 1 

_^j> 

_d 


w _ P 4> 03 ^ ^ ^ »v 

„ Sm'o'H.tg 3j = " d< 2 o e o 

l&^S 2 “ E J * H-5 3 

*? 2 “** "^3 

c n % c * _ 2 £ - o S o o 

s ?«! s "2 -* H *5 e 


S 'S bp 
O U «S 


cd S 

3 « O 


Lrj w cr --i _ o > ~ O 03 ^ w 

S|«1 3-3- sf? J I g 

J! K ° -* 7 | = u 1"3 2 i'rt 

— + - 3 « R ° « — a, c 

s-sgfc-ssaossSo 

fl^Q-ocoJuOo) o o 0 


c 

C 


= -; : 
bC P ’" 
g £ 

J s c 

s "o M 

W U u- 

^ 'o 0 

2 - * v 

I go 

ou-3 

S o 5 

; <tf 

1 .£ ^ * r 
! ^ er 

> 2 03 

1 "3 £ v 

; > d g 

; - -T 3 -2 

!.G ii o 

| K a3 &. 

, 5 *w S 

* g o •- 
1 » 

is S-J 

! h jg J 

* -*> t> 

; o t3 2 

. M C 03 

: tJ nj (j 


•2 ^ S 

C w as 

3 _ 03 

O «J g 

I- M M 

w 3 “3 

5 ° 'S 

d 2 © 

« xS c 

03 «J P 

3* c o) 
$ « J3 


W3 ed u " 

— •> H - 

C r3 L os 

cj o ^ a 
j- S m >-> 

O n ^ , 
a 03 w - 

Ed S *>' 

.S u M 

*5 -a 

S-Ji” 

^ "S bo 

£ £ 6 E 

E ,g 

*> £ £ « ; 
>-« t- ftj 

* ff.Sl' 

z C 2 e 

T3 •-« ft, 

r* Z1 03 

5 to -P JS 

^ .2 c t — 

J5 — O 

«E h ’S 5 

-"'5*1 J 

n u Zr -P 

03 c » 03 

■O «1 e-o 

O ■** !• , 

£ o 05 -i 

- §• 5> s 

-«s u 2. 

S ° g s 

-2 03 <d ft, 

S -o c JS 

y fc-t >» CJ 

'-0 ° JS 

.5 to £ ^ 

CO *S 03 ? 

„ S - no 

3 2 « g 
-2 w j- - 
P « V -t? 

H -t? ? — * 


O £ *7? f *-■ 

u M3 -5 , -5 £ -, 

v '3 s e '* n 

^3 m P > *“0 < 

- 81 "j « » * i 

° S H E-I " 2 J 

c .s? £ . -o -s - 

g s S 2 10 o ' 

So „ s u 

O. J3 .— o 3 ® o 

-*-» r~i >-• O ^ 

»> . • «> CJ • « 

a — <= Z e * 

•o= M 2 % "« ' 
c II “ «j w: w “o ’ 
a5 ^ E -2 03 O 

. SICV * 60 -D g 

nO'^Q C - c 

03 • P ^ 03 ® 

JS 03 7? M «j > 

JS JS ° « O rl 

a» -** > js tr 1 

H S - S - u . ■ 

O- > * o 2 £*>' ■ 

o ■ 5 - c ^ _ 

1211 fi^: 

oj *■ ■> o *2 P 
> 2 T3 Jx W 5 . 

■ 5 ja p. *5 *i 5 

° -. <fl d “ ^ 

s » o 5 2^. 

■O JK SL 5 5, t S 

_ O s s -o ° ' 

£ S 2 w £ « o ' 

2 .- 3 c scr- 

2 1 s = u- ° 8 

2 2 E « o “ o 

-• a ^ x *o 

“ > 2 = C — 

^ S « 2 -2 > « 

. u _c 3 aS oS 

i o O* "5 a c > 


E c *^3 -a 

C3 0) JX 03 

*» s T3 E 

- P X O -- 

« | i 6 & 

j 5 w u i> 

T3 *P ^ O *- 
C « 

ft, 03 gj 

J £ as - 

3 u- tO c «— 

8 S7*a ° 

2 ^ 12 
o) r\ <- 
& tfi 03 o 

o o u kO . 

*p -* ° flS 

c CO 43 Vi 3 
oS 03 W5 O O 

e s 2 12 -O 
■S* u ~ "2 ® 

M a^s .is .2 

l 03 O 7» 3 

! - O > c rTS 

I 3 7 ^ J 

1 g 1? ■** 3 

. ^ 03 fi « ^ 

^* S I 1 *2 * 2 -d 

: 2 ^ S « « S 
i £^ s>»£ 13 

5 8 - - S-a O 

i »J.S J» s 

= K " £. It .2 


> £-5 s 

« S. . - 3 
s ?:«PJ g 
bj >a &o 5 o 



Table S. Maximum mkgniludcs of weighted aerodynamic terms. 
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